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Abstract 

Q-conditional symmetries (nonclassical symmetries) for a general class of two- 
component reaction-diffusion systems with constant diffusivities are studied. Using the 
recently introduced notion of Q-conditional symmetries of the first type (R. Cherniha J. 
Phys. A: Math. Theor., 2010. vol. 43., 405207), an exhaustive list of reaction-diffusion 
systems admitting such symmetry is derived. The form-preserving transformations for this 
class of systems are constructed and it is shown that this list contains only non-equivalent 
systems. The obtained symmetries permit to reduce the reaction-diffusion systems under 
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of the form 

03 



l/"~) ■ study to two-dimensional systems of ordinary differential equations and to find exact so- 

lutions. As a non-trivial example, multiparameter families of exact solutions are explicitly 
constructed for two nonlinear reaction-diffusion systems. A possible interpretation to a 

^ ■ biologically motivated model is presented. 
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1 Introduction 

The paper is devoted to the investigation of the two-component react ion- diffusion (RD) systems 
form 

u t = d 1 u xx + F(u,v), ,_ 

Vt = d 2 v xx + G(u,v). 

where u = u(t,x) and v = v(t,x) are two unknown functions representing the densities of 
populations (cells), the concentrations of chemicals, the pressures in thin films, etc. F and G 
are the given smooth functions describing interaction between them and environment, d\ and d 2 
are diffusivities assumed to be positive constants. The subscripts t and x denote differentiation 
with respect to these variables. The class of RD systems (CQ) generalizes many well-known 
nonlinear second-order models and is used to describe various processes in physics, biology, 
chemistry and ecology (see, e.g., the well-known books [H [21 El IH [5] ) . 

Nevertheless the search for Lie symmetries of the class of RD systems (pQ) was initiated 
about 30 years ago [6], this problem was completely solved only during the last decade because 
of its complexity. Now one can claim that all possible Lie symmetries of (CQ) were completely 
described in [ZHEHS]. 



The time is therefore ripe for a complete description of non-Lie symmetries for the class 
of the RD systems ([JJ). However, it seems to be extremely difficult task because, firstly, sev- 
eral definitions of non-Lie symmetries have been introduced (nonclassical symmetry [TJ ITU] , 
conditional symmetry [TT| [12] . generalized conditional symmetry [131 EH] etc.), secondly, the 
complete description of non-Lie symmetries needs to solve the corresponding system of deter- 
mining equations, which is non-linear and can fully be solved only in exceptional cases. 

Hereafter we use the most common notion among non-Lie symmetries, non-classical symme- 
try, which we continuously call the Q-conditional symmetry following the well-known book [11] 
and our previous papers [151 US]- It is well-known that the notion of Q-conditional symmetry 
plays an important role in investigation of the nonlinear evolution equations because, having 
such symmetries in the explicit form, one may construct new exact solutions, which are not 
obtainable by the classical Lie machinery. However, for a complete description of such sym- 
metries, one needs to solve the corresponding non-linear system of determining equations that 
usually is very difficult task. Thus, to solve the Q- conditional symmetry classification problem 
for the class of RD systems (JTJ), one should look for new constructive approaches helping to 
solve the relevant nonlinear system of determining equations. A possible approach was recently 
proposed in [T7] and is used in this paper. 

It can be noted that the diffusion coefficient d\ in system (TjQ) can be omitted without losing 
of generality because the simple substitution 

t -»■ t/dx, F -> -diC 1 , G -> -d 2 C 2 

reduces the system to the form 

u xx = u t + C 1 (u,v), .. 

v X x = dv t + C 2 (u,v), 

where d = 4 3 -. Thus, we consider system (J5J) in what follows. 

The paper is organized as follows. In section 2, two different definitions of Q-conditional 
invariance for the class of RD systems ($Z§ are presented and the system of determining equations 
is derived. The theorem giving the complete description of Q-conditional symmetries of the first 
type is proved. In section 3, the form-preserving transformations for the class of RD systems 
02]) are constructed and applied to the RD systems derived in section 2. In section 4, the Q- 
conditional symmetry obtained for reducing of the RD systems to the ODE systems are applied. 
Examples of finding exact solutions are presented together with a possible interpretation for 
population dynamics. Finally, we summarize and discuss the results obtained in the last section. 

2 Conditional symmetries of the RD systems 

Here we use the definition of Q-conditional symmetry of the first type for the RD systems (see 
|17] for details). It is well-known that to find Lie invariance operators, one needs to consider 
system (J2J) as the manifold A4 = {S\ = 0, S 2 = 0} where 

Si = u xx -u t -C x {u,v), 

S 2 = v xx - dv t - C 2 (u, v), 



in the prolonged space of the variables: t, x, u, v, u t , v t ,u x , v x , w a 



u xt ,v xt ,u tt ,v tt . According 



to the definition, system (|2]) is invariant under the transformations generated by the infinitesimal 
operator 



Q = £°{t,x,u,v)d t + £, 1 (t,x,u,v)d x + r] 1 (t,x,u,v)d u + r] 2 (t,x,u,v)d v 
if the following invariance conditions are satisfied: 



(3) 



QSi = Q (u xx -u t - C x {u, v)) 

2 2 

QS 2 = Q (v xx - dv t - C 2 (u, v)) 
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The operator Q is the second prolongation of the operator Q, i.e. 
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where the coefficients p and o with relevant subscripts are expressed via the functions £°, £*, rj 1 
and rf by well-known formulae (see, e.g., [HI [181 US])- 

Hereafter the listed above differential operators act on functions and differential expressions 
in a natural way, particularly Q (u) = (,°u t + ^u x — rj 1 and Q (v) = £°v t + ^v x — if. 

Definition l.pZ] Operator ([3]) is called the Q-conditional symmetry of the first type for the 
RD system (J2J) if the following invariance conditions are satisfied: 



QSi = Q (u xx -u t - C x (u, v)) 



0. 



Mi 



QS 2 

2 



Q (v xx 

2 



dv t — C 2 (u,v)) 



0. 



Mi 



where the manifold .Mi is either {Si = 0, S 2 = 0, Q (u) = 0} or {Si = 0, S 2 = 0, Q (v) = 0}. 

Definition 2. Operator (|3j) is called the Q-conditional symmetry of the second type, i.e., the 
standard non-classical symmetry for the RD system (|2j) if the following invariance conditions 
are satisfied: 



QS X = Q (u xx -u t - C x {u, v)) 

2 2 ' 

QS 2 = Q (v xx - dv t - C 2 (u, v)) 



M-2 



Mi 



0. 



0. 



where the manifold M 2 = {Si = 0, S 2 = 0, Q (u) = 0, Q (v) = 0}. 

Remark 1. It is easily seen that M. 2 C M.i C Ai, hence, each Lie symmetry is automatically 
the Q-conditional symmetry of the first and second type, while each Q-conditional symmetry 
of the first type is one of the second type (non-classical symmetry). 

Remark 2. To the best of our knowledge, there are not many paper devoted to search of 
Q-conditional symmetries for the systems of PDEs [20l I2TI [22| [23l 121] . One may easily check 
that Definition 2 was only used in all these papers. 



Statement. Let us assume that 

X={hi(t,x)v + ho{t,x))d v (4) 

(hereafter hi(t,x) and ho(t,x) are the given functions) is the Lie symmetry operator of the 
RD system (EJ) while Q\ is the known Q-conditional symmetry of the first type, which was 
found using the manifold Mi = {Si = 0, S2 = 0, Q (u) = 0}. Then any linear combination 
C\Q\ + C2X (C\ and C2 7^ are arbitrary constants) produces new Q-conditional symmetry 
of the first type. 

Application of definition 2 for finding Q-conditional symmetry (non-classical symmetry) 
operators of the RD system (J2J) leads to a complicated system of determining equations ( DEs) 
(see system (19) in [17j), which seems to be extremely difficult for exact solving. 

It turns out that application of definition 1 leads to essentially simpler system of DEs, 
which can be fully integrated. Here we present the result under the restrictions (° / and 
d 7^ 1 (the cases £° = and d — 1 must be investigated separately). Thus, the system of DEs 
corresponding to the manifold Mi = {Si = 0, S 2 = 0, Q (u) = 0} takes the form 

e x = e u = e v = e u = e v = o, (5) 

11 2 2 2 n / c \ 

Vv = Vuu = Vuu = Vvv = Vuv = °> ( 6 ) 

2« u + (^-l)« = 0, (7) 

Hu + £ = 0, (8) 

27& + det=0, (9) 

a£ - &° = 0, (10) 

v'cl + V 2 Cl + (2£ - vDC 1 = V l - V l (11) 

^ + ^ + (2^^^^ + (l-4^ + ,L-^. (12) 

Note that there is no any need to solve the similar system of DEs corresponding to the 
manifold Ai\ = {Si = 0, S2 = 0, Q (v) = 0} because the discrete transformations u — > v, v — ¥ u 
transform each symmetry found using .Mi to one corresponding to the manifold Ai*. 

It should be also noted that we find purely conditional symmetry operators, i.e., exclude 
all such operators, which are equivalent to Lie symmetry operators described in [TJEj. Having 
this aim, we use the system DEs for search Lie symmetry operators (see [H] for details): 

(13) 

(14) 

(15) 

2^ + ^ = 0, (16) 

Ho + d£ = °> ( 17 ) 

V'CI + r, 2 Cl + (2£ - vDC 1 = V l x - rjl (18) 

V'CI + V 2 C 2 V + (2£ - vl)C 2 = vl ~ dvl (19) 





■£ = 


St; 




e 1 = 


0, 


vl = 


-vl- 


~ 'luu 


Ivv 


= 0, 




2^- 


-$ 


= 0, 









Comparing DEs (15|)- (IT2"|) with (IT5|) -( ITU|) one concludes that r] 2 ^ is the necessary and sufficient 
condition, which guarantees this property. 

Now we need to solve the nonlinear system (l5l)- ([T2l . Obviously equations (jSJ) and (JSJ) can 
be easily integrated: 

r/ 1 = r x (t, x)u + p l (t, x), rf = q(t,x)u + r 2 (t,x)v + p 2 {t,x), 

where £°(£), £ x (t, x), q(t,x), r k (t,x), p k (t,x) (k = 1,2) are to-be-determined functions. Thus, 
substituting (120]) into fl7|)-( TT2|) . one obtains the nonlinear system of PDEs: 

2£V + £V- l)g = 0, (21) 

K + e t = 0, (22) 

2r* + d& = 0, (23) 

2£ - tf = 0, (24) 

( r i u + pi)C7i + (gu + r ^ + v *)Cl + (2£* - r 1 )^ 1 = (r^ - rf )u + p^ - P*, (25) 

( r i u + tf)Cl + (g M + r 2 v + p 2 )C 2 + (2£ - r 2 )C 2 = qC 1 + ^±^g(l - rf) + (26) 



0"L - *1> + (9a» - d*)« + pL - rfp : 



* • 



to find the functions £°(£), £*(£, a;), q(t,x) ^ 0, r k (t,x), p k (t,x). In other words, all possible 
Q-conditional symmetries of the first type are easily constructed provided the general solution 
of system (T2TT) — ff26|) is known. 



Theorem 1 The nonlinear RD system (dp with d ^ 1 is invariant under the Q-conditional 
operator of the first type (GJ) if and only if one and the corresponding operator have the forms 
listed in Table 1. Any other RD system admitting such Q-conditional operator is reduced to one 
of those from Table 1 by the local transformations 

t-*C 1 t + C 2 , 

x — > C%x + C4, , , 

u -> C 5 e C6t u + C 7 , [ ' 

v -)- C 8 e Cf)t v + C10, 

with correctly- specified constants Ci, I — 1, . . . , 10 and/or by adding a Lie symmetry operator of 
the form ([J]) 

Sketch of proof. To prove the theorem one needs to solve the nonlinear PDE system (J2"T|) - 
(|26|) with restriction q(t, x) 7^ 0. We remind the reader that C l and C 2 should be treated 
as unknown functions. As follows from the preliminary analysis (see equations (j25p and (|2"6"1) 
involving the functions C l and C 2 ), we should examine 6 cases: 

(1) r 1 = r 2 = p 1 = 0, 

(2) r 1 = r 2 = 0, p 1 ^ 0, 

(3) r 1 = p 1 = 0, r 2 ± 0, 



(4) r 1 = 0, pV 0, r 2 ^ 0, 

(5) r 2 = 0, r V 0, 

(6) r 1 ^ 0, r 2 ^ 0. 

Solving system f l2T]) - fl26l) in each case one obtains the list of Q-conditional symmetries of the 
first type together with the correctly-specified functions C 1 and C 2 . Note that the symmetry 
operators have the different structures depending on the case. 

Let us consider case (1) in details. Equations (F25]) and (T2"6"]) take the form 

(qu + p 2 )C 2 + 2HC 2 = qC 1 + (q xx - dq t )u + p 2 xx - dp 2 . [ ^> 

Differentiating the first equation of (F2SJ) with respect to x, one arrives at the equation 
(q x u + p 2 )Cl = 0, which lead to the requirement C* = 0. In fact, if q x ^ then immediately 
C\ = 0. If q x = then equation (|2T|) produces £* = 0, hence, C* = 0. Thus, the first equation of 
system (1281) takes the form £, x C l = and two subcases £j; ^ and £* = should be examined. 

The general solution of (]2"g]) with £* ^ is 

C 1 = 0, C 2 = exp ( - -^- 2 v)g(u) + ^--^ + ^-^, (29) 

where g(u) is an arbitrary (at the moment) function. Because the function C 2 doesn't depend on 
t and x, equation ( 129]) with g{u) ^ immediately produces the restrictions q = «i£*, p 2 = «2^i, 
where a.\ and a?2 are arbitrary constants. Differentiating equation ( )24l) with respect to x, one 
obtains £, xx = 0. So q x = ai^L = 0, however, this contradicts to the assumption £* 7^ 0. The 
remaining possibility g(u) = leads to the linear RD system (T5]). 

Now we examine the subcase £* = 0, i.e., £* = Ai = const. The general solution of ( 128]) 
takes the form 

^1 w n ^2 g/M + faxx - dq t )u + p 2 xx - dp 2 

C = /(«), C = q~^f v + g{u), (30) 

where f(u) and g(u) are arbitrary (at the moment) functions. 

If f(u) is an arbitrary function then we obtain p 2 — f3q (/3 — const ) Hence C 2 = j^rlv + av + 

g(u), where a = qxx ~ qt . Having this, we use renaming ^2!l — > f{u) and solve the overdetermined 
system 

qxx-dqt _ 

2q*+\ 1 (d-l)q = 0. 

Thus, the system of DEs (121]) - (126]) is completely solved (under above listed restrictions !) and 
we obtain the conditional symmetry operator 

„ o xo /Ai(l-d) \ 2 (l-d) 2 -4a \. 

Q = d t + \ x d x + A 2 exp [-^ ! -x + -± -j tj (u + f3)d v , 



where Ai and A2 7^ are arbitrary constants, of the RD system 

u xx = u t + (u + (3)f(u), 

v X x = dv t + f(u)v + av + g(u). 

Finally, using the simple transformation 
u —¥ u — ft, 



(31) 



(32) 



one sees that it is exactly case 6 of Table 1. 

To complete the examination of case (1) we look for the correctly-specified function f(u), 
which satisfies (1301) without the restriction p 2 = (3q. Indeed, if one finds the differential con- 
sequences of the second order (see equation for C 2 ) then C 2 X = 0, C 2 t = and two algebraic 
equation to find the function f(u) are obtained: 

(q t p 2 - qp 2 )f = ((q xx - dq t )u + p 2 xx - dp 2 )(q t u + p 2 )- 

((q X x - dq t ) t u + (p 2 xx - dp 2 ) t )(qu + p 2 ), 
(q x p 2 - qp 2 x )f = {{q xx - dq t )u + p 2 xx - dp 2 )(q x u + p\)- 

{{q X x ~ dq t ) x u + (p 2 xx - dp 2 ) x )(qu + p 2 ), 

Thus, f{u) = a\ + a^u + a^u 2 provided p 2 7^ f3q. Substituting this expression into (!30|) and 
making the standard routine, one arrives at case 9 of Table 1 if 0:3 7^ and case 10 if «3 = 0. 

Cases (2)-(6) were treated in the similar way and the results are listed in Table 1. It should 
be noted that several local transformations ( (1321 is the simplest example) were used to reduce 
the number of cases and simplify structures of the relevant RD systems. These transformations 
can be presented in the general form (T2TJ) . 

The sketch of proof is now completed. ■ 

Table 1. Q-conditional symmetry operators of the RD system <^/ with d ^ 1. 





C\u,v) 


C 2 {u,v) 


Q 


1. 


u/(w) 


UI = u (v — u) 


dt + aud u + a((l — k)u + kv)d v , 


2. 


uf(uj) 


u{g{ui) + ail — d) In u+ 
/(w)lnn), w = uex.p(-f i ) 


dt + aud u + a(u + v)d v , a^O 


3. 


u/(w) 


g{uo)+u{f{uo)+a{l-d)) 
uj = uexp(u — v) 


dt + aud u + ot(u + l)d v , a^O 


4. 


/(«) 


e u g(uj) - f(u) - a(l - d) 

ui = e~ u {u + v) 


dt + ad u + a{u + v — l)d v , a/0 


5. 


/(«) 


uf(u)+g(u) + (l-d)u 

ui = u 2 — 2v 


dt + d u + ud v 



6. 


uf(u) 


vf(u) + g(u) + av 


d t + ^d x + qud v , A2/O 
q = A 2 exp f Aix + ^-^0 


7. 


/(«) 


(u + v)( g(u) + a ln(n + v) J — 


9 t + \exp(-%t)(u + v)d v , 

A/0 


8. 


/(«) 


a(w + u) — /(u) 


<9 t + A(u + t>)ft„ A/0 


9. 


«1 + «2"U+ 

u 2 


<7(u) + ut> 


ft + f^ + fow + P 2 )^. 

q = (pi (i) exp (Aix) , (px ^ 

P 2 = ((^1 + a 2 )(pi - dipij exp (Aix) 


10. 


«1 + OL 2 U 


3(«) + a 3 w 


ft + T^ftz + M + ^ft, 

g = A 2 exp ^Aix + - J — 3 tj , 

#L = dpf + a 3 p 2 - aiq, A 2 + 


11. 


«4 ln(u + u) 


a 3 -u + (a 3 - a 2 )n- 
a4ln(u + v) 


a t + (V.(x)exp(^^t)-a)a- 
ft) + a "(\- 2 d) (~» + u )ft' a 1 a 2 a 4: ^ 


12. 


ct^it + v 


a 3i ,+ i(l + d)^+ 
aiulnu + «4ti 


9* + <p 2 (t)(ud u + uft,) - tp 2 (t)ud v , 

U 2 (f 2 (t) / 


13. 


a 2 U 


(X3V + 04ti + ti fc 


<9t + \ 2 ud u + (993 (t)u + X 2 kv)d v , \ 2 <P3 ¥" 


14. 


aiulnu 


a 3 i> + aivlnn+ 

1 
a2^ d 


ft + T^dx + \ 2 e~ ait ud u + (qu + f e"" 1 ^)^, 

9 = A 3 exp (Aix + ^t + ^A 2 e- Q i*) 

CK1A2A3 7^ 


15. 


aiulnu 


a 3 v + «iw In u + a^u 


ft + A 2 e- Ql 'ud u + (¥> 4 (*)« + (if e- ait + 
A 3 )u)<9„, aiA 2 ^4 / 


16. 


a±ulnu 


(X3V + aif lnu+ 

1 

CK4U + a2ti d 


ft + A 2 e- ai 'n5 u + (tp 5 (t)u + ^e-^V)ft, 
01020^2^5 / 


17. 


aiulnu 


031; + a\dv lnu+ 
«i (1 — d)u In u — a^u 


ft + A 2 e- ait ua u + (am + {X 2 e- ait + ai)v)ft, 
aiA 2 7^ 



18. 


Oil + 2 U 


03^ + Inn 


dt + jz^d x + X 2 ud u + (qu + p 2 )^, 

g = A3 exp ^Aix + -^ ^ i '-tj , 

p 2 xx = dp 2 + CX3P 2 - a±q + A 2 , A 2 A 3 / 


19. 


CKl + Ct 2 U 


013V + Inn + a±u 


<9 t + \2ud u + ((fe(t)u +p 2 )d v , a±X 2 tpQ 7^ 0, 
pL = d Pt + «3P 2 - ai^e(i) + A 2 


20. 


a 2 u 


013V + ulnu 


d t + A 2 u<9„ + (<P7(t)u + \2v)d v , X 2 cp 7 ¥" 


21. 


OL\ + a 2 U+ 

u 2 


uv + as 


d t + X(v + (p 8 (t)u + a 2 ip 8 (t) - dip 8 (t))d v , 


22. 


OL\ + « 2 U 


03V + u 2 


dt + j^d x + P 1 d u + {qu + p 2 )d v , 

q = <£> 9 (i)exp(Aix), V3g 7^ 

PL = ^4 + «3P 2 + (d - l)qp l - aiq, 

p 1 = |((Af + a 2 - a 3 )c^9 - #9) exp (Aix) 


23. 


ai 


0:3V + e u 


dt + jz^d x + A 2 3 M + (qu + A 2 w + p 2 )d v , 

vlx = d Pt + a 3P 2 - (A 2 (l - d) + ai)q, 

g = A 3 exp(Aix+^^t), A3 + 


24. 


OL\ 


0:3V + 014U+ 
e u 


3* + A 2 <9 U + ((pio{t)u + A 2 v + p 2 )<%, a 4 A 2 99io / 0, 
p 2 ^ = dp 2 + a 3 p 2 - (A 2 (l - d) + ai)q + a 4 A 2 


25. 





a\(u + v)x 
\n(u + v) + a 2 


dt + T 2 i(d u - d v ) + Aexp(-^t)(u + «)$,, 
«ia 2 A 7^ 


26. 





u 2 


2tft + x9, + ^ exp(-|J)9 u + (^ exp(-g)u+ 
2v+p 2 )d v , d = 3, ?4 = 3p 2 + ^ exp(-f^), 


27. 





u 5 


t 2 9 t + ta<9 x - ^ud u + (Aexp(-g)u- 
faa + a v)ft , d-5, A/0 



In Table 1, the functions t/>(x), y?i(£), tp 2 (t), <P4(t), <P5(t), <P8(t), and v?9(0 are the general 
solutions of the equations 



a\a 2 



v a 4 (l — a) y 

Af + a 2 A 2 + «i 



dV?i - (2X 1 + a 2 )v?i + 



d 



<Pi = 0, 



d</3 2 - («2 - a 3 + (1 - d)ip 2 )ip 2 - a ± ip 2 = 0, 



d^ + (a 3 + Hd - l)e- Qlt )^ 4 - a 4 (A 3 + A 2 ^i- Ql< ) = 0, 
# 5 + (« 3 + X 2 (d - l)e- ai4 )^ 5 - a 4 A 2 ^i- ai4 = 0, 

«1 a 3 

dy3 8 - a 2 y? 8 + -r^s + -» = °> 



and 



d</3 9 - (Ai(l + d) + a 2 (l - d) - "3)^9 + (A? - a 3 Ai - a 2 (a 2 - "3)^9 = 5 
respectively. The functions 



fs{t) 



<Pio(t) 



A 3 exp ^-0,+Wi-o q + ^f^ , if « 2 ^ « 3 - A 2 (l - d), 

-^^t + A 3 , i/ a 2 = a 3 -A 2 (l-d); 



d 



A 3 exp(-^t) + ^, ifa 3 ^0, 



IT "> 



s^at, if a 3 = 0. 



Finally, the function y?6(£) = y? 3 (£) at A; = 0, while y?7(£) = y3 3 (£) at A; = and a 4 = 1. Hereafter 
the upper dot index denotes differentiation with respect to the variable t. 

3 Form-preserving transformations of the RD systems 

A natural question is: Can we claim that 27 systems listed in Table 1 are inequivalent up to 
any local substitutions (not only of the form (1271) !)? It turns out that the answer is positive. 
To present the rigorous proof of this, we used the notion of the set of form-preserving point 
transformations introduced in [25] and now extensively used for Lie symmetry classification 
problems (see, e.g. [26| [27]). Note that finding these transformations for systems of PDEs is 
a difficult problem because of technical problems occurring in computations and there is no 
many results for systems (paper [28] is one of the first presenting an explicit result for a class 
of systems) . 

The form-preserving transformations present the most general and correctly-specified form 
of local substitutions, which can map some equations from a given class to other those belonging 
to the same class. They contain as particular cases the well-known equivalence transformations 
and discrete transformations, which maps each equation from the class to another one from 
this class, used in the well-known Ovsiannikov method of Lie symmetry classification. Here we 
construct such transformations with the aim to show that Table 1 cannot be shortened. 

Theorem 2 An arbitrary RD system of the form (TJ|) with d / 1 can be reduced to another 
system of the same form 

w yy = w T + F 1 (w,z), ,_ 

Zyy = \Z T + F 2 (W,Z), 

10 



by the local non-degenerate transformation 

r = a(t, x,u,v), y = b(t, x, u, v), 
w = <p(t,x,u,v), z = ip(t, x,u,v ), 



(34) 
(35) 



if and only if the smooth functions a, b, ip and ip are: 

a = a(t), b = (3{t)x + j(t), 

ip = f(t) exp ( - ^(/3x 2 + 2 7 x)) w + P(t, x), ( 36 ) 

i/j = g{t) exp ( - ±0x 2 + 2^?x))v + Q(t, x), 

where the functions a(t), f3(t), f{t), g(£),7(£), P(t,x), and Q(t,x) are such that the equalities 

a = (3 2 , X = d, (37) 

fi 2 F 1 ((p, ip) = VuC l {u, v) + Lf xx -ft- 2—ip xu , (38) 



ft - 2 ip xu , 

fu 



(5 2 F\f^) = ip v C 2 ( Uj v) + ip xx - dipt ~ 2^ 



(39) 



take 'place provided /3fg ^ 0. 



Proof. First of the all we note that each non-degenerate transformation ( 1341) - ( 1351) should 
satisfy the conditions 



Ai 



a x a t 
b x bt 



7^0, 



A, 



fu fv 
Ipu tpv 



7^0, 



(40) 



which are used to prove the theorem. 

Let us choose an arbitrary RD system of the form ([2]) . The main idea of the proof is based 
on substituting the expressions for u xx , v xx , u t , v t using the formulae (134"|) and (135]) into this 
system and on analysis conditions when the system obtained is equivalent to system (1331 . The 
expressions for the first-order derivatives have the form 



tp x - a x w T - b x w y a v w T + b v w y - cp l 

lyx Q'xZt O x Zy a v z T ~\~ O v Zy ip v 



u, 



U t 



a u w T + b u w y — ip u a v w T + b v w y — (p. 

O'uZt T" O u Zy %p u CL v Z r T" O v Zy 1JJ V 

(ft - a t w T - b t w y a v w T + b v w y - (p v 

ipt - d t z T - b t z y a v z T + b v z y - ip v 



a u w T + b u w y - ip u a v w T + b v w y - <p v 

UuZf T" U Zy lp U Qj V Z t "T OyZy tfjy 



The expressions for the second-order derivatives are very cumbersome, however, it can be noted 
that they contain the derivative w TT and w ry . Because r is a new time- variable we conclude 
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that the coefficient next to w TT and w Ty must vanish otherwise system ( 13 3 p are not obtainable. 
These coefficients vanish if and only if the equalities take place: 

a x = a u = a v = b u = b v = =>- a = a(t), b = b{t, x). (41) 

Moreover, taking into account (I4"D"1) . the restriction a b x ^ is also obtained. 

Having the set of equalities (14J1 . the expressions for u xx and u t can be essentially simplified, 
namely: 

"xx ^ 2 W ra Aa ^ y?/ -| ^ Wy -£2 6yT 

{■>Px'fiv-^vVx)x^2~{^-2)x{lpxfv-1pvVx) (42) 

A| > 

M * = lk {^{awr + 6 t u><, - (ft) - (p v (az T + b t z y - ip t )) . 

Substituting (142]) into the first equation of ([2]). Omitting the full expression of the equation 
obtained, we note that one contains the terms 

r-\ ¥vb 2 x , a , .... i) v b 2 x/ a , 

0) - -^\Zyy ~ pZrh W -^-( w w _ ^)' 

while other terms don't depend on z yy ,Wyy,z T , and w T . 

Now there is two possibilities. If the first equation of ([2]) is transformed into the first one 
of (133 p then we immediately obtain 

a = b 2 x , y v = 0. (43) 

If the first equation of ([2]) is transformed into the second one then the conditions 

a = Xb 2 x , ip v = (44) 



must be satisfied. It turns out that conditions ff44|) lead to the result, which is obtainable from 
(136]) by the discrete transformations u — > v and v — > u. 

Let us consider conditions ( 143]) . Taking into account ( 140]) and (p v = 0, the restriction 
ipuipv j^ springs up. 

On the other hand, b(t,x) = (3{t)x + j(t), follows from ( 143]) where /3 and 7 are arbitrary 
smooth functions. Thus, the first equation from (137]) is derived. 

Substituting (143]) into expressions for u xx and u t (see formulae ( 142]) ). one obtains 

„, _ /3 2 „,, o P<p xu I 2ifi x tp xu -ip u (p xx <fi uu ( O \2 

U xx ~ ^Wyy Z^Wy+ ^2 ^3 lPW y ip x ) , ^ 

u t = ^(aw T + 0x + j)wy-tp t ). 

Since the first equation of system (1331) doesn't contain the terms w y and w y , we should vanish 
the relevant coefficient, namely: 

2§£xu _|_ /3x+7 _ Q 

W« = 0. 
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The general solution of this system can be easily constructed so that obtains 

if = fit) exp ( - -^0x 2 + 27x)) M + P(t, x), (46) 

where f{t) ^ and P{t, x) are arbitrary functions at the moment. Thus, the first, second 
and third equations from fl36l) are derived. Moreover, substituting fj4"5]) and (J4"6]) into the first 
equation of system ([2]), we arrive at the equation 

w yy = w T + -^{C (u,v)-- — ). (47) 



Now one realizes that (jl"7)) coincides with the first equation of system (IH3T) iff condition 
takes place. 

The analogous routine involving the second equation of system ([2]) leads to the condition 
a = ^/3 2 =>• A = d (see (1431 ). the function if) of the form (136]) and equation (|39|) . 

The proof is now completed. ■ 

Consequence 1. The set of transformations ( 12"T|) arising in theorem 1 is a subset of 
form-preserving transformations (I3"6"j) . 

Consequence 2. If the nonlinear RD system of the form (JT|) is transformed to another 
one from this class, say, to the system 

w*. =d\u* x , x ,+F*{u\v*), 
v;* = d* 2 v** x * + G*(u*,v*) 

by a local substitution then they have the proportional diffusivities. Moreover, there are two 
linear combinations for the reaction terms F and F*, and for G and G* resulting ot\U + a 2 and 
a^v + 0:4, respectively (here a^, A; = 1, . . . , 4 are correctly-specified constants). 

Roughly speaking, consequence 2 says that the locally-equivalent RD systems have the same 
structure up to additive terms a\u + a 2 and a^v + 0J4. At the first sight, there are some systems 
in Table 1 satisfying this consequence, for example in cases 18 and 19. However, according to 
consequence 2, the term a 3 u arising in the second equation of the RD system (see case 19) 
cannot be removed by any local substitution. We have carefully checked all cases listed in 
Table 1 and concluded that there aren't any locally-equivalent systems therein. 

Thus, we have shown that the list of RD systems presented in Table 1 cannot be reduced 
(shortened) by any local substitution. 

4 New exact solutions and their possible interpretation 

It is well-known that using the known Q-conditional symmetry (non-classical symmetry), one 
reduces the given system of PDEs to a system of ODEs via the same procedure as for classical 
Lie symmetries. Since each Q-conditional symmetry of the first type is automatically one of the 
second type, i.e., non-classical symmetry, we apply this procedure for finding exact solutions. 
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Thus, to construct an ansatz corresponding to the given operator Q, the system of the linear 
first-order PDEs 

Q (u) = 0, Q (v) = (48) 

should be solved. Substituting the ansatz obtained into the RD system with correctly-specified 
coefficients, one obtains the reduced system of ODEs. 

Let us construct exact solutions of the non-linear RD system listed in the case 1 of Table 1, 
when the system and the corresponding symmetry operator have the form 

u t = u xx -uf{u), , 4g , 

dv t = v xx - u k g(u) - u(f(u) + q(1 - d)), u = u~ k {y - u) 

and 

Q = d t + aud u + a((l — k)u + kv)d v , (50) 

In this case system ( I48p takes the form 

Ut = a n * 4. fc ( 51 ) 

v t = a(l — k)u + akv 

and its general solution produces the ansatz (the functions u and v depend on two variables t 
and x): 



u = ip(x)e at 

v = ?p(x)e kat + tp(x)e 



kat i , n f„\„cet V*2) 



where <p(x) and ip(x) are new unknown functions. Substituting ansatz fl52|) into (jUJ), one 
obtains so called reduced system of ODEs 



if;" = tp k g(uj) + akdip, u = ipip 



k\(, ,\ i ^,Lj„/, , , _ j„„-k (^3) 



Because system fl53l) is non-linear (excepting, of course, some special cases) it can be inte- 
grated only for the correctly-specified functions / and g. We specify / and g in a such way, 
when the RD system in question will be still non-linear (otherwise the result will be rather 
trivial). Thus, setting f(ui) = 70;^ —a, g(u>) = (3u, ( (3 and 7 are arbitrary non-zero constants, 
the RD system takes the form 

Ut = u xx - j(v - u)k + au, , , 

dv t = v xx — 7(t> — u)~k — (3v + ((3 + ad)u, 

while the corresponding reduced system is 

y> = (p + akd)i). K ' 
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The general solution of (I55"|) can be easily constructed: 

^ x) = y f(f^ {x)dx ) dx + C3X + C4 , (56) 



{C\ exp(/ix) + c 2 exp(— /ix), if /i 2 = /3 + afcrf > 0, 
cicos(ra) + c 2 sin(ra), if i/ 2 = — (f3 + afcd) > 0, (57) 

C\X + c 2 , if /3 + a/cd = 0. 

Thus, substituting (J5S|) and ( 157)) into (152")) . the 4-parameter family of solutions for the non-linear 
RD system (1541 is constructed. 

Hereafter we highlight the solutions satisfying the zero Neumann boundary conditions, which 
widely arise in biologically motivated boundary- value problems. Hence, setting c 3 = c 4 = 0, k — 
|, if) = c\ cos(ux), one obtains the solution 



u = — 7<^ (cos 2 (ra) + 6) cos(ra)e at , 
v — c\ cos(z/x)es ai + it. 



It can be noted that this solution satisfies the zero Neumann boundary conditions 

u x \ x =o = 0, v x \ x=0 = 0, u x \ x=j e. = 0,v x \ x= jz = 

on the interval [0, j-], where j G N. 

Let us set f(u) = — (oi + 6a;), g(o;) = (a(l — d) — oi)a; (hereafter a, a\ and 6 are arbitrary 
non-zero constants) in (|49p . hence, it takes the form 



Mt = M„ + OlW — 6« 2 k + &WM 1 k 

dv t = v xx — a 2 v — bu 2 ~ k + bvu 



^ ™ I £i/i t n i -L rC V 



where a 2 = a(l — d) — a\. The corresponding reduced system of ODEs is 

V " + tyip 1 -" + (a x - a)<p = 0, 

<//' = (afcd + a 2 )V- l J 



Nevertheless we have not constructed the general solution of system (1591) . its particular solution 
was found by setting if) = —5, 5 ^ 0. In this case, the first-order ODE 



• = ±^(a - a^ + ^- + c„ tt = jj^l^, * + 2 (60) 

for the function (p is obtained ( the value k = 2 is special and leads to ODE <p' = 
±a/(« - a{)(p 2 + 26(51nv? + ci). 

If ci 7^ then the general solution of ( 160)) can be expressed via hypergeometric functions. 
Here we present the solution for (160)) with ci = 0, k ^ 0: 



fl( t an*(*^(ar±c*))+l)) ", a 1>a , 

ip(x) = { ) \-i ( 61 ) 

^(tanh^^p^xic,))-!)) , oi<a, 
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Figure 1: Exact solution ( 1631) with k = 0.5, 5 = 6, d = 4, ai = 5, 6 = 3. 



( ^ = 2bs )' w hi cri seems to be the most interesting. Note that another arbitrary constant 
can be removed by the trivial substitution x ± C2 — > x. 

Now we rewrite system (155)) setting t> — >■ — 1> with the aim to obtain a biologically motivated 
model. So the system takes form 



Ut = u xx + u(ai — bu 1 k ) — bvu 1 k , 
dv t = v xx + v(—a 2 + 6n 1_fc ) + bu 2 ~ k , 



(62) 



where all coefficients (excepting k) should be positive. ( 162)) can be treated as a prey-predator 
model for the population dynamics. In fact, the species u is prey and described by the first 
equation. Its population decreases proportionally to the predator density v. The natural birth- 
dead rule for the prey is u(a\ — bu 1 ^ k ) and can be treated as a generalization of the standard 
logistic rule u(a% — bu) (see, e.g., [2]). The similar arguments are also valid for the second 
equation. The model should involve also the zero Neumann boundary conditions (zero-flux on 
the boundaries), which indicate that both species cannot widespread over the globe but occupy 
a bounded domain. 

Using ( 1521) with v — > —v, ip — —& and ( IBT)) with oi > a we construct the exact solution 



u 



,5 (tan 



2( kyJa-L-a 



x) + r 



.at 



Se 



akt 



/3(tan 



2( ky/ai — a 



x) + l) 



(63) 



at 



of ( 162)) . It turns out that solution can describe interaction between prey and predator on the 
space interval [0,/], (here / = k J* 3 _ n , j £ N) provided 



fci/ai- 



d' ~ V 26 



a: 



Oi 



d(ife -1) + 1 



< 1. 



(64) 
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One easily checks that solution ( 16 3 j) is non- negative, bounded in the domain Q = {(t, x) G 
(0, +oo) x (0, 1)} and satisfy the given zero Neumann boundary conditions, i.e. 

u x \x=o — 0, v x \ x= q = 0, u x \ x= i = 0, v x \ x= i = 0. 



As example we present this solution ( 16311 with the parameters satisfying the restrictions ( 1641) 
in Fig. [TJ This solution can describe such type of the interaction between the species u and v 
when both of them eventually die, i.e. (u, v) — > (0, 0) if t — > +oo. 

5 Conclusions 

In this paper, Q-conditional symmetries for the class of RD systems (T5]) (that is equivalent to 
the class of systems ([I])) and their application for finding exact solutions are studied. Following 
the recent paper [17], the notion of Q- conditional symmetry of the first type was used for these 
purposes. The main result is presented in theorem 1 giving the exhaustive list of RD systems of 
the form ([I]) with d x ^ d 2 (the case di = d 2 should be analyzed separately), which admit such 
symmetry. It turns out that there are exactly 27 locally-inequivalent RD systems admitting 
the Q-conditional symmetry operators of the first type of the form ([3]) with ^° ^ (the case 
£° = should be analyzed separately). To show local non-equivalence of the systems listed 
in Table 1, we proved theorem 2 describing the set of form-preserving point transformations 
for the class of RD systems ([2]). Note that all the operators found are inequivalent to the Lie 
symmetry operators presented in [7J E] because the necessary and sufficient condition, which 
guarantees this property, was used. 

The Q-conditional operator listed in case 1 of Table 1 was used to construct the non-Lie 
ansatz and to reduce two nonlinear RD systems to the corresponding ODE systems. Solving 
these ODE systems, the two-parameter families of exact solutions were explicitly constructed for 
the RD systems in question. Moreover, application of the exact solutions for solving the prey- 
predator system ([62]) was presented. It turns out that the relevant boundary value problem 
with the zero Neumann conditions can be exactly solved and the solution can describe the 
densities of two interacting species. 

The work is in progress to construct conditional symmetries for multicomponentKD systems. 
In particular case, a wide list of the Q-conditional symmetries of the first type for the three- 
component diffusive Lotka-Volterra system is presented in [29] . 

Finally, we point out that this paper is a natural continuation of the recent paper [16J, 
where RD systems with non-constant diffusivities were examined. 
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